An extended-phase-space dynamics for the generalized nonextensive thermostatistics 



O 

o 



o 



J. S. Andrade Jr., M. P. Almeida, A. A. Moreira, and G. A. Farias 

Departamento de Fisica, Universidade Federal do Ceard, 

Campus do Pici, Caixa Postal 6030, 60455-760 Fortaleza, Ceard, Brazil 

(February 1, 2008) 

We apply a variant of the Nose-Hoover thermostat to derive the Hamiltonian of a nonexten- 
sive system that is compatible with the canonical ensemble of the generalized thermostatistics of 
Tsallis. This microdynamical approach provides a deterministic connection between the generalized 
nonextensive entropy and power law behavior. For the case of a simple one-dimensional harmonic 
oscillator, we confirm by numerical simulation of the dynamics that the distribution of energy H 
follows precisely the canonical g-statistics for different values of the parameter q. The approach is 
further tested for classical many-particle systems by means of molecular dynamics simulations. The 
results indicate that the intrinsic nonlinear features of the nonextensive formalism are capable to 
generate energy fluctuations that obey anomalous probability laws. For q < 1 a broad distribution 
of energy is observed, while for q > 1 the resulting distribution is confined to a compact support. 
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Since the pioneering work of Tsallis in 1988 fm where 
a nonextensive generalization of the Boltzmann-Gibbs 
(BG) formalism for statistical mechanics has been pro- 
posed, intensive research has been dedicated to de- 
velop the conceptual framework behind this new thermo- 
dynamical approach and to apply it to realistic physical 
systems. In order to justify the generalization of Tsal- 
lis, it has been frequently argued that the BG statistical 
mechanics has a domain of applicability restricted to sys- 
tems with short-range interactions and non (multi)fractal 
boundary conditions g. Moreover, it has been recalled 
that anomalies displayed by mesoscopic dissipative sys- 
tems and strongly non-Markovian processes represent 
clear evidence of the departure from BG thermostatistics. 
These types of arguments have been duly reinforced by 
recent convincing examples of physical systems that are 
far better described in terms of the generalized formal- 
ism than in the usual context of the BG thermodynamics 
(see || and references therein). It thus became evident 
that the intrinsic nonlinear features present in the Tsal- 
lis formalism that lead naturally to power laws represent 
powerful ingredients for the description of complex sys- 
tems. 

In the majority of studies dealing with the thermo- 
statistics of Tsallis, the starting point is the expression 
for the generalized entropy S q , 



1 - / [f{x)]Hx 
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where k is a positive constant, q a parameter and / is the 
probability distribution. Under a different framework, 
some interesting studies [|4| have shown that the param- 
eter q can be somehow linked to the system sensibility 
on initial conditions. Few works have been committed to 
substantiate the form of entropy (fil) in physical systems 
based entirely on first principles pLpfl. For example, it 
has been demonstrated that it is possible to develop dy- 
namical thermostat schemes which are compatible with 



the generalized canonical ensemble m . In a recent study 
by one of us |8J, a derivation of the generalized canoni- 
cal distribution is presented from first principle statistical 
mechanics. As a consequence, it is shown that the partic- 
ular features of a macroscopic subunit of the canonical 
system, namely, the heat bath, determines the nonex- 
tensive signature of its thermostatistics and therefore its 
power law behavior. More precisely, it is exactly demon- 
strated in M that if one specifies the capacity of the heat 
bath as 
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where q is a constant, 1/(3 oc kT, and T is the temper- 
ature, the generalized canonical distribution maximizing 
(0) for the Hamiltonian H of the system is recovered, 



f(H)<x[l + P{q-l){E-H)] — 



(3) 



Here, E is a conserved quantity and denotes the energy of 
the extended system (system-)- heat bath). Equation (g) 
provides a very simple but meaningful connection be- 
tween the generalized g-statistics and the thermodynam- 
ics of nonextensive systems. It is analogous to state that, 
if the condition of an infinite heat bath capacity is vio- 
lated, the resulting canonical distribution can no longer 
be of the exponential form and therefore should not fol- 
low the traditional BG thermostatistics. In the present 
letter, we will show how the conjecture proposed in M 
can be used to develop a variant of the original Nose- 
Hoover thermostat p 10 that is consistent with the q- 
thermostatistics. We will then validate the technique by 
applying it to the cases of a simple harmonic oscillator 
and a classical many-particle system. 

We consider a system of N particles having coordi- 
nates x[, masses rrii and potential energy $(a/). As in 
the extended system method originally proposed by Nose, 



here we also introduce an additional degree of freedom 
through a variable s which will play the role of an exter- 
nal heat bath, acting as to keep the average of the kinetic 
energy at a constant value. In practice, this is achieved 
by simply rescaling the real variables in terms of a new 
set of virtual variables 



where, 
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where A is a rescaling exponent, and (x' i ,p' i ,t') and 
(xi,pi,i) are the real and virtual coordinates, momenta 
and time, respectively. At this point, we postulate that 
a generalized Hamiltonian for the extended system can 
be written as 



H q (x,p,p B ,s)=^2 
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where the first two terms on the right side represents the 
energy of the physical system that is free to fluctuate 
JUj. The virtual variable p s also has a real counterpart, 
Ps = Ps/s x , and has been introduced to allow for a dy- 
namical description of the variable s. More precisely, the 
third term p 2 /2Q corresponds to the kinetic energy of 
the heat bath and the parameter Q is an inertial fac- 
tor associated with the motion of the variable s. The 
last term of the Hamiltonian (0) is a power law poten- 
tial in s. As we show next, it provides the essential link 
between the concept of extended-phase-space dynamics 
and the generalized canonical ensemble. 

We start by considering the quasi-ergodic hypothesis 
and writing the time average of a given quantity A(x', p') 
as 



~* = \ 



dx dp dp s ds —5(H q — E) (6) 



with Z= ill I dyLdpdp s dss- l 8(H q -E) , 

where Z is the analogous to a microcanonical normaliza- 
tion factor for the generalized Hamiltonian (g) . The fac- 
tor s _1 in (0) accounts for the variability in the real time 
sampling during the extended-phase-space dynamics || . 
Transforming the virtual momenta p and coordinates x 
back to real variables, changing the order of integration 
and rewriting the volume element as dxdp = s gX dx'dp' 7 
where g is the number of degrees of freedom, we obtain 



A = j f fdx'dp'A (I d Ps ds s s x - 1 5{H q 
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If we now make use of the property of the 8 function, 
i5[/j(s)] = S(s — s )/h'(s ), where s is the zero of h, it 
follows that 
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H = H(x',p') = '£?i- + $(x>). 
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By integration with respect to p s we get 



Z\2 1 J v 2' 7 



x / dx! dp' A[l + or/(E - H)] 
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where B is the beta function, B(z,w) = J Q t z : (1 — 
tY'^dt. Finally, if we define a = f3(q + l)/2g\ and 
7 = 2g\(q — l)/(q + 1), the generalized canonical aver- 
age is recovered 
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A = — j j dx'dp'A[l+0(q-l)(E-H)]T^ (11) 
dx'dp' [l+/3(q-l)(E-H)]T^) , 



Z 
with Z = 



and we have thus proved that, under conservation of the 
extended Hamiltonian Eq. (^), the fluctuations in the 
energy H(x' , p') of the physical system should be consis- 
tent with the canonical formulation of the nonextensive 
q-thermostatistics. 

It is possible to confirm the validity of this approach 
with a simple realization of the generalized thermostat 
scheme. We consider an extended system composed of 
a single one-dimensional harmonic oscillator coupled to 
a heat bath whose thermal capacity obeys essentially 
Eq. (0). From Eq. (pi), such a system can be described 
by the following extended Hamiltonian: 



H q (x,p,p s ,s) 
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Here m — 1 for simplicity and we choose to set A = 2 
because the nonlinear dynamics for this case when q = 1 
(i.e., for the BG thermostatistics) has been shown to be 
sufficiently chaotic to generate average properties of the 
canonical ensemble [fL2l. From (O) and the scaling re- 
lations (11), we obtain the equations of motion for the 
extended system in the real phase space, 
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A fifth order Runge-Kutta subroutine is then used to 
numerically solve this set of nonlinear differential equa- 
tions. To ensure the conservation of energy H and the 
stability of integration, all runs have been performed with 
10 8 time steps of At' = 10~ 4 each. The density maps 
shown in Figs. 1(a) and (b) for q = 0.8 and 1.2, respec- 
tively, provide clear evidence that the dynamics of both 
systems fills space. For q > 1 [see Fig. 1(b)], the ac- 
cessible phase-space lies in a compact set, whereas the 
phase-space support for q < 1 [see Fig. 1(a)] is infinite. 
The former situation is compatible with the necessary 
cut-off condition on energy for q > 1 j|. In Fig. 2 
we show the logarithmic plot of the distributions of the 
transformed variable \ — 1 + P(q ~ 1)(-® ~~ H), where 
H = p' 2 /2+x' 2 /2, for three different values of the param- 
eter 7 = 4(g— l)/(g+l) corresponding to q — 0.7, 0.8 and 
0.9. Indeed, we observe in all cases that the fluctuations 
in x follow very closely the prescribed power law behav- 
ior, p(x) ck X : an d therefore confirm the validity 
of our dynamical approach to the generalized canonical 
ensemble. As shown in Fig. 3, the simulations performed 
for q > 1 are also compatible with the expected scaling 
behavior. However, instead of the long-range tail ob- 
tained for the case q < 1, a rather unusual power law 
with positive exponent is observed. 

Now we focus on a more complex application of the 
thermostat scheme introduced here. The basic idea is to 
simulate, through molecular dynamics (MD), the nonex- 
tensive behavior of a classical many-particle system. For 
completeness, we start by rewriting the expression for 
the extended Hamiltonian (H) in terms of the usual q- 
thermostatistics parameters (i.e., q and (3) 
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where ln q (s) = (s^ 1 - l)/(q - 1) @. From Eq. |[J), it 
is then possible to derive the equations of motion for any 
value of q and any type of effective potential of interac- 
tion. We consider a cell containing 108 identical parti- 
cles that interact through the Lennard-Jones potential, 
$(Axij) = 4e[(cr/Axy) 12 - (a/Axij) 6 ], where Axy is 
the distance between particles i and j, e is the minimum 
energy and a the zero of the potential. Distance, energy 
and time are measured in units of a, e and (ma) 2 /e, re- 
spectively, and the equations of motion are numerically 
integrated using a predictor-corrector algorithm [pL3| . 

Compared to the previous example of a single har- 
monic oscillator, the complexity of the many-particle sys- 
tem hinders a quantitative prediction of the statistical 
behavior of its energy fluctuations. Because the exact 
form or even a plausible approximation of the density of 
states tt(H) is difficult to obtain in this case, we restrict 
ourselves to the qualitative analysis of the resulting en- 
ergy distribution p(H) oc £l(H)f(H). Furthermore, we 
performed additional simulation tests with different num- 
ber of particles and physical conditions to confirm that 



the MD system is always led to unstable trajectories in 
phase-space whenever the value of q is set to be smaller 
than a given threshold (? m i n . In spite of these limitations, 
however, the results shown in Fig. 4 clearly indicate the 
tendency for a broader distribution of energy when q < 1 
(we set q to be slightly larger than q m [ n sw 0.9940 in this 
case). For q > 1, on the other hand, the resulting dis- 
tribution of energy is notably more confined than the 
Maxwellian distribution obtained for q — 1 (see Fig. 4). 

In summary, we have shown that the essential features 
of the generalized canonical distribution can be captured 
with a proper extension of the standard Nose-Hoover 
thermostat. To the best of our knowledge, this is the 
first time that a Hamiltonian approach to the nonex- 
tensive q-thermostatistics leads explicitly to the obser- 
vation of a power law behavior (q < 1). We thus be- 
lieve that the microdynamical formalism presented in this 
work can provide a deterministic link between the gen- 
eralized entropy conjecture Eq. (|lj) and the concept of 
Levy flights R,H. Finally, the methodology introduced 
here is flexible enough to accommodate the description 
of other nonextensive systems of physical significance. 
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FIG. 1. (a) Density plot of the harmonic oscillator dy- 
namics subjected to the generalized thermostat scheme for 
q = 0.8. The initial conditions are [x'(0) = 0.5, p'(0) = 0.5, 
s(0) = 1.0, p' s (0) = 0.0] and the thermostat parameters have 
been set to a = 1.0 and Q — 1.0. (b) Same as (a) but for 
q = 1.2. 
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FIG. 2. Logarithmic plot of the distributions of the trans- 
formed variable \ for q — 0.7 (circles), 0.8 (squares) and 0.9 
(triangles). From right to left, the three straight lines with 
slopes —3.33, —5.0 and —10.0 correspond to the expected 
power law behavior p(x) oc X ■ 




FIG. 3. Logarithmic plot of the distributions of the trans- 
formed variable \ for q — 1.1 (circles), 1.2 (squares) and 1.3 
(triangles). From right to left, the three straight lines with 
slopes 10.0, 5.0 and 3.33 correspond to the expected power 
law behavior p(x) oc x ■ 
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FIG. 4. Logarithmic plot of the energy distributions for 
q = 0.9941 (circles), 1.0 (full circles) and 1.1 (triangles). In 
all three cases, the MD simulations have been performed with 
108 particles, f3 — 0.2, and a density of 0.1 particles/cr 3 . 



